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Using fundamental-measure density functional theory we investigate entropic wetting in an asym- 
metric binary mixture of hard spheres with positive non-additivity. We consider a general planar 
hard wall, where preferential adsorption is induced by a difference in closest approach of the dif- 
ferent species and the wall. Close to bulk fluid-fluid coexistence the phase rich in the minority 
component adsorbs either through a series of first-order layering transitions, where an increasing 
number of liquid layers adsorbs sequentially, or via a critical wetting transition, where a thick fllm 
grows continuously. 

PACS numbers: 68.08.Bc, 64.70.Ja, 82.70.Dd 



Studying the interfacial properties of liquid mixtures is 
of significant fundamental and technological relevance [1] . 
Bulk liquid-liquid phase separation, which can arise at or 
close to room temperature, is usually associated with rich 
phenomenology of interfacial behaviour at a substrate. 
Gaining a systematic understanding of how the differ- 
ent types of intermolecular and of substrate-molecule in- 
teractions induce phenomena such as wetting, layering, 
and drying at substrates constitutes a major theoreti- 
cal challenge. Relevant for surface adsorption of liquids 
are Coulombic and dispersion forces, but also solvent- 
mediated and depletion interactions which occur in com- 
plex liquids. Arguably the most important source for 
the emergence of structure in dense liquids is the short- 
ranged repulsion between the constituent particles; this 
may stem from the overlap of the outer electron shells 
in molecular systems or from screened charges or steric 
stabilization in colloidal dispersions. 

Hard sphere fluids form invaluable reference models for 
investigating the behaviour of liquids at substrates. Both 
the pure |11 [3] and binary [4] hard sphere fluids are rel- 
evant, the latter playing an important role when adding 
e.g. electrostatic interactions in order to study wetting 
of ionic liquids at a substrate [5]. The most general bi- 
nary mixture is characterized by independent hard core 
distances between all different pairs of species, and is 
referred to the non-additive hard sphere (NAHS) model. 
Here the cross species interaction distance can be smaller 
or larger than the arithmetic mean of the like-species di- 
ameters. The NAHS model gives a simplified representa- 
tion of more realistic pair potentials, i.e. charge renormal- 
isation effects in ionic mixtures in an explicit solvent in- 
duce non-additive effective interactions between the ions 
[S] . It is also a reference model to which attractive or re- 
pulsive tails can be added [7]. The Asakura-Oosawa-Vrij 
(AOV) model of colloids and non-adsorbing polymers |8] 
is a special case where one of the diameters (that of the 
polymers) vanishes. 

It is surprising that the wetting behaviour of the gen- 
eral NAHS model is largely unknown, given the funda- 
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FIG. 1: (a) Illustration of the asymmetric NAHS model with 
positive non-additivity. The solid boundaries represent the 
hard cores of the small and big species. The dotted line repre- 
sents the non-additive hard core between unlike species, which 
here is attributed only to smaller particles, (b) Three exam- 
ples of general planar hard walls. The additive wall treats the 
two species equally, while the 6-type and s-type walls have 
properties similar to the big and small particles, respectively. 



mental status of the model. In this Letter we address this 
problem and consider the NAHS fluid at a general, non- 
additive hard wall. We find a rich phenomenology of in- 
terfacial phase transition, including two distinct types of 
surface transitions: one is layering, where the adsorption 
of one of the phases occurs through a number of abrupt 
jumps, and the other is critical wetting, where the thick- 
ness of the adsorbed film grows continuously when vary- 
ing the statepoint along the bulk fiuid-fluid binodal. Via 
changing the wall properties a crossover between these 
transitions occurs. 

The binary NAHS model is defined by the pair poten- 
tials Vij{r) = oo for r < aij and otherwise, where i,j 
= s,b refers to the small and big species, respectively, 
ass and abb are the diameters of the small and big parti- 
cles, respectively, and r is the center-to-center distance. 
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The cross-species diameter is asb = 5(1 + A)(i7ss + (^bb), 
where A > — 1 measures the degree of non-additivity, 
see Fig. [ija) for an illustration of the length scales. The 
model is characterised by the size ratio, q — OssI Obb ^ 1, 
and by A. In this Letter we restrict ourselves to the 
asymmetric size ratio q = 0.5, and to positive non- 
additivity A = 0.2, as a representative case. We re- 
late A to a length scale via d = ^(ctss + o'f,b)A = 
asb — ■^{o'ss + crbb), where here d — O.Sass- The statepoint 
is characterised by two partial bulk packing fractions, 
iji — naf^pi/G, where pi is the number density of species 
i. We define a general planar hard wall via the external 
potentials Ui(z) = 00 if z < k, and otherwise; here z is 
the distance between the wall and the particle center, and 
li is the minimal distance of approach of species i = s,b. 
Clearly the origin in z is irrelevant, so the only further 
control parameter is the wall offset, 61 = lb — Is- For addi- 
tive hard sphere mixtures it is common to set li — au/i; 
for our model parameters this results in 61 — 0.5ass- Be- 
sides this 'additive wall', two further special cases are 
shown in Fig.jlJ^b). The 6- type wall has properties similar 
to the big particles so that it sees these with their 'intrin- 
sic' size lb = Cbb/'i, but sees the small particles with their 
'non- additive' size Is = crss/'2 + d, such that 61 = 0.2(7^3. 
We expect that the bigger particles adsorb more strongly 
to this wall. Conversely, the s-type wall has properties 
similar to the small particles, so that it sees these with 
their 'intrinsic' size Is — Ossj'^-, and sees the big parti- 
cles with their 'non-additive' size lb = 'Jbbl'^ + <i-, so that 
61 = 0.8(7ss- Thus, one expects the small particles to 
adsorb more strongly. 

We investigate the inhomogeneous NAHS fluid using 
a fundamental measure density functional theory [SJ [TU] . 
Comparison of theoretical results to Monte Carlo simula- 
tion data for bulk fluid- fluid phase diagrams [HI [TU] , par- 
tial radial distribution functions [HI HJ and density pro- 
files in planar slits [T^] indicates very good quantitative 
agreement. We obtain equilibrium density distributions 
Piiz) from the grand potential functional, by 
numerical solution of 6^l/6pi{z) = 0, i = s,b. To calcu- 
late coexisting (bulk or surface) states we use the equal- 
ity of the chemical potentials ps, IJ-b-, and in the two 
phases. The NAHS functional [3] features both a large 
number of terms and a large number of convolutions that 
take account of the non-locality. Therefore the accurate 
calculation of density profiles close to phase coexistence, 
and close to interfacial transitions, is a challenging task. 

For q = 0.5 and A = 0.2 the DFT predicts fluid- 
fluid phase separation with a critical point at 77^ = 0.049, 
rjb = 0.151 - see Fig. 2(b). We start with the b-type wall, 
which we find does indeed preferentially adsorb the big- 
ger particles. For 6-rich statepoints the preferred species 
is already at the wall and no surface transitions occur. 
For s-rich statepoints at bulk coexistence, but far from 
the bulk critical point, we find that the small particles 
dominate the region close to the wall, but that there is 
a small amount of adsorption of the bigger particles. To 
illustrate this, see the pair of density profiles, Ps{z) and 
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FIG. 2: (a) Pairs of density profiles, Ps{z) (left panel) and 
Pb{z) (right panel), of the NAHS fluid with q — 0.5 and 
A — 0.2 at a 6-type wall and at bulk coexistence on the s- 
rich side of the phase diagram. The shaded regions represent 
the range of profiles possessing n = or 1,2,3,4 or 5 ad- 
sorbed 6-rich layers, and the region where the adsorbed film 
becomes infinitely thick. The solid lines represent specific ex- 
amples from the middle of each range. The inset shows the 
adsorption of each species, Fi, as a function of rjs- (b) The 
corresponding phase diagram in the {risjrjb) plane. There is 
a series of layering transitions that intersect the bulk binodal 
(*) and descend into the s-rich one phase region, ending in a 
surface critical point (o). For clarity only the first two tran- 
sitions are shown in full, while the remaining transitions are 
represented only by their intersection with the bulk binodal. 
The inset shows the location of the first layering and the 'wet- 
ting' transitions in relation to the bulk critical point (•). 



Pb{z), furthest from the bulk critical point in Fig. [2ja). 
Reducing rjs along the binodal in the direction towards 
the bulk critical point, there occurs a series of discon- 
tinuous jumps of the density profiles. The first jump 
corresponds to the big particles displacing the small par- 
ticles from the wall and forming a layer at a distance abb 
away from the wall, see Fig.[2](a). Each subsequent jump 
corresponds to the adsorption of an extra 6-rich liquid 
layer at the wall. Using the coexistence criteria we have 
located five distinct layering transitions. Beyond the fifth 
transition we find that the layer rich in the big particles 
becomes macroscopically thick. We discuss the possible 
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FIG. 3: Same as Fig.|2|a), but for the s-type wall. The coex- 
istence curve is traced on the 6-rich side of the phase diagram. 
As the wetting critical point is approached the smaller parti- 
cles strongly adsorb at the wall, replacing the bigger particles 
and growing a thick film. Below the wetting critical point 
the film is infinitely thick. The inset shows the adsorptions, 
Ti , against the difference in the packing fraction of the small 
species from its value at the wetting critical point, rj* — rjs, 
on a logarithmic scale, where i]g — 0.0043. 



nature of this transition below. The inset to Fig. [2]^ a) 
shows the adsorption, = J dz [pi{z) — pi{oo)], of each 
species i = s,b as a, function of 775. Each plateau rep- 
resents the range of statepoints along the binodal which 
have a particular number of adsorbed layers. The forma- 
tion of the infinitely thick layer corresponds to Ff, jump- 
ing to -l-oo, and F^ to —00. The layering transitions are 
first-order surface phase transitions, characterised by a 
range of coexisting states. In Fig. [2j[b) we plot the co- 
existence lines of the first two transitions in the {r]s,r]b) 
plane. We find that the layering transitions intersect the 
bulk binodal [T^ and that they lie very close to the bin- 
odal on the s-rich side of the phase diagram. Each tran- 
sition terminates at a surface critical point, where the 
jump in Ti vanishes. The first layering transition, where 
the big particles strongly adsorb at the wall and form the 
first layer, is the largest both in terms of the change in 
the adsorptions and its size on the phase diagram. Each 
subsequent transition is smaller than the previous one. 

We next turn to the s-type wall. As this preferentially 
adsorbs the smaller particles, tracing the bulk coexistence 
curve on its &-rich side is interesting. For statepoints far 
from the bulk critical point, we find that there is some 
adsorption of the small particles, but that big particles 
dominate the region close to the wall, see the pair of 
density profiles furthest from the bulk critical point in 
Fig. 3, where Pb{z) exhibits oscillatory decay that indi- 
cates high-density packing effects. Increasing rjg along 
the binodal in the direction of the bulk critical point, we 
find that the small particles start to adsorb more strongly 
at the wall, replacing the big particles. On moving fur- 
ther towards the bulk critical point, a thick film rich in 
the small particles grows. No jumps are observed and the 



FIG. 4: Value of rjs at (i) the intercept of the layering transi- 
tions with the bulk binodal (solid lines) and (ii) the location 
of the critical wetting transition critical point (dashed line), 
as a function of the scaled wall offset 51 /ass- As SI /ass is 
increased from 0.2 (&-type wall) the layering transitions move 
along the binodal towards the bulk critical point, located at 
•q"'^ ~ 0.05. At SI /ass — 0.27 the layering transitions coa- 
lesce and the surface transition becomes critical wetting. As 
51 increases towards the additive case, the critical wetting 
transition approaches the bulk critical point. Increasing 51 
further, the wetting critical point moves to the other side of 
the binodal. The inset shows the s-type wall wetting critical 
point (a) in relation to the bulk critical point (•). 



thickness increases continuously (and reversibly) with the 
state point up to a wetting critical point, beyond which 
the film is infinitely thick, see Fig. [3j Hence we con- 
clude that this wetting transition is critical. In such a 
case the adsorption can be shown [13j [14] to diverge as 
Ti cx log(|?7* — 77s|) on the mean-field level, where 77* is 
the value of rjs at the wetting critical point. We find the 
value of ry* by fitting F^ to its asymptotic form. The inset 
to Fig. [3] compares the adsorptions to the asymptotic log- 
arithmic form. The location of the wetting critical point, 
?7* = 0.0043 is shown in relation to the bulk binodal in 
the inset to Fig. |4j 

We next vary the wall offset parameter, SI, between the 
two cases discussed above. Starting with the 6-type wall, 
SI /ass = 0.2, and increasing SI we find that the location 
of the layering transitions moves towards the bulk critical 
point. In Fig. |4] we show the value of r]s at each of the 
intersections of a layering transition and the bulk bin- 
odal as a function of SI. The jump in adsorption at each 
layering transition becomes smaller and the extent of the 
line in the phase diagram becomes shorter (not shown). 
Decreasing SI further, we find that at SI /ass — 0.27 the 
individual layering transitions bunch up and become in- 
distinguishable from each other. For smaller SI there is 
a single continuous wetting transition, where the thick- 
ness of the adsorbed 6-rich layer grows logarithmically, 
in a similar manner to the behaviour at the s-type wall 
described above. We establish the location of the surface 
critical point by fitting F.; to its asymptotic form and plot 
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the value of rj* at the wetting critical point in Fig. [4] In- 
creasing SI further results in the location of the wetting 
critical point moving further along the bulk binodal to- 
wards the bulk critical point so that at Sl/a-gs — 0.43 the 
wetting transition critical point coincides with the bulk 
critical point, and the wall is neutral such that neither 
species is preferentially adsorbed at the wall. As SI is in- 
creased beyond 0.43 we find that the wetting transition 
moves to the 6-rich side of the phase diagram. The addi- 
tive wall, SI/ ass = 0.5, has a critical wetting transition, 
but located very close to the bulk critical point. As SI is 
increased, the wetting critical point moves further along 
the bulk binodal, so that we return back to the s-type 
wall, Sl/ass — 0.8. 

In order to ascertain the generality of our findings, we 
have investigated the trends upon changing the model 
parameters. For size ratio q — 0.5 and vanishing wall 
offset, SI = 0, we find layering transitions far from the 
bulk critical point for a range of non-additivity param- 
eters A = 0.1,0.2,0.5. Adjusting SI towards the case of 
the additive wall, the layering transitions move towards 
the bulk critical point. We also investigated symmetric 
mixtures with q = 1 and A = 0.1. Clearly, for the ad- 
ditive wall, (5 = 0, there is no preferential adsorption at 
the wall and hence no layering transitions. Introducing 
preferential adsorption via a non- vanishing wall offset, 
SI = 0.1,0.2,0.3, layering transitions occur, and these 
move away from the bulk critical point upon increas- 



ing 51. 

In summary, we have shown that the NAHS model ex- 
hibits both layering and critical wetting transitions de- 
pending on the hard wall offset parameter. We expect 
this wetting scenario to be general and to occur in a large 
variety of systems where steric exclusion is relevant. A 
set of layering transitions had been previously found in 
the AOV model at a hard wall [TS]. In these studies the 
wall parameter is equivalent to the 6-type wall. As in 
these previous papers the existence of an infinite number 
of layering transitions is a possibility within our mean- 
field DFT treatment. The effects of fluctuations would be 
to smear out the higher-order layering transitions to pro- 
duce a final 'wetting' transition as found here. A change 
from a first-order to a critical wetting transition is not 
uncommon [14 . What is remarkable here is that tri- 
critical behaviour can be induced in a purely entropic 
system by merely changing a non-additive wall parame- 
ter, SI. Moreover, the NAHS model is much less special 
than the AO model, as here both species (not only the 
AO colloids) display short-ranged repulsion and hence 
packing effects. In future work, it would be interesting 
to see the effects of non-additivity on wetting in charged 
systems where first-order and critical wetting transitions 
occur [5]. 
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